Interfacial Slip in Sheared Polymer Blends 
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We have developed a dynamic self-consistent field theory, without any adjustable parameters, 
for unentangled polymer blends under shear. Our model accounts for the interaction between 
polymers, and enables one to compute the evolution of the local rheology, microstructure and the 
conformations of the polymer chains under shear self-consistently. We use this model to study the 
interfacial dynamics in sheared polymer blends and make a quantitative comparison between this 
model and Molecular Dynamics simulations. We find good agreement between the two methods. 
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Self-consistent field (SCF) theory has been very useful 
for modeling and understanding the equilibrium prop- 
erties of inhomogeneous polymer fluids. In lattice SCF 
theories, Markovian statistics for chains represented by 
lattice random walks are coupled with a mean field ap- 
proximation for the free energy to determine variation 
of equilibrium properties across interfaces 0. In con- 
trast, interfacial dynamics and rheology of inhomoge- 
neous polymer fluids driven out of equilibrium by applied 
stresses is still not well understood. Since experimental 
evidence indicates that the interfacial properties are criti- 
cal to the evolution of the morphology of polymer blends, 
understanding the interplay between the interfacial dy- 
namics and the polymer conformation at the interface is 
of great fundamental importance, and efficient compu- 
tational models of these systems are essential for many 
industrial and biological applications. 

We have recently developed a novel dynamic self- 
consistent field (DSCF) theory by combining the lat- 
tice random walk formalism of Scheutjens and Fleer's 
static SCF model and a convective-diffusive lattice 
gas model for simple compressible fluids ■ As in static 
SCF theory, the distribution of free segments interact- 
ing with a self-consistent field representing adjacent con- 
nected segments and walls is approximated by a product 
of free segment (one-body) probabilities. The effect of 
flow on ideal chain conformations is modeled with FENE- 
P dumbbells, and related to stepping probabilities in a 
random walk. Free segment and stepping probabilities 
generate statistical weights for chain conformations in a 
self-consistent field, and determine local volume fractions 
of connected segments. Flux balance across a unit lattice 
cell yields mean-field transport equations for free segment 
probabilities and momentum densities. Diffusive and vis- 
cous contributions to the fluxes arise from segmental hops 
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modeled as a Markov process. Hopping transition rates 
depend on the changes in the free energy, reflecting seg- 
mental interactions, kinetic energy, and entropic contri- 
butions accounting for chain deformation under flow. It 
is the first dynamic mean-field theory for inhomogeneous 
polymer fluids that couples self-consistently the effects of 
flow, segmental and wall interactions, and compressibil- 
ity on chain conformations, composition, transport and 
rheology in interfacial regions at the Kuhn length scale, 
which is necessary to resolve interfacial properties of typi- 
cal polymer interfaces. This makes it different from other 
recently proposed approaches Q • In this Letter, we focus 
on the interfacial dynamics of a symmetric polymer blend 
in the unentangled regime. We use the DSCF model to 
predict the velocity slip, chain conformations, and the 
reduction of shear viscosity across the interface between 
the coexisting polymer phases, and compare our results 
with molecular dynamics (MD) simulations. 

We consider a blend of two immiscible homopolymer 
species A and B consisting of linear chains of N and N B 
segments, respectively. The fluid is confined between two 
parallel walls normal to the z-axis in a Cartesian coor- 
dinate system at temperature T. The interface between 
the two polymers is parallel to, and is centered midway 
between the two walls. The walls move at equal and op- 
posite velocities along the x-axis. We discretize the space 
between the two walls into triangular lattice layers par- 
allel to the walls, stacked in a FCC lattice arrangement. 
The fluid system considered here is invariant under lattice 
translations in the xy-plane, and we simplify the DSCF 
equations utilizing this symmetry (see |J for details). 

Let Pf, (/>£ and gf be the free segment probabil- 
ity, segmental volume fraction and segmental momen- 
tum density of species a (ev=A or B) in each layer i 
(i=l, 2, • ■ • , L), respectively. The free segment probabil- 
ities and the segmental momentum density g;=g^ + g,^ 
evolve according to 
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which are obtained by applying mass and momentum 
conservation laws to a rectangular control volume of 
w=a 3 /v / 2 centered at a site in layer i. Here Cfi±i * s the 
locally averaged segmental friction coefficient of species 



a (see for details) and Ui=giW / (m A cpf + m B (j)f), 
where m a is the segmental mass of species a, is the 
mass averaged velocity at a site in layer i. The tensor 
—ef in Eq. represents the elastic contribution to the 
stresss and is defined in 0. Let D? l+1 =k B T/N a (^ +1 
and vf i+1 —C t f i+1 a 2 /2^m a be the locally averaged self- 
diffusion coefficient and kinematic viscosity of species a, 
and 4> be the overall average polymer volume fraction, the 
diffusive free segment probability current and the viscous 
stress tensor at the mid-plane between layers i and i + 1 
in Eqs. (JTJ and © are given by 
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In Eqs. J2J and ijljl. (p is the Kawasaki transition rate 
function <p(x)=2/(l + e x ) which satisfies local detailed 
balance, and A{H")=( K H^, 1 ) — {H") is the free energy 
change due to a hop of a segment of species a from layer 
i to i + 1. The local contribution to the free energy of a 
segment of species a in layer i is 



(i 



"a/3/ 



X^B«^f»-X?(*il+*i£) 



/3=A,B 



In det 



i>fk B T 
2N a 



Tr I - 



3 gig 



(5) 



where xab is the segment-segment interaction parameter 
for species A and B, and x" is t ne segment- wall inter- 
action parameter of a segment of species a. In Eq. J5J, 
the double angular brackets represent summation over all 
nearest neighbors of a site in layer i and Sf is the sec- 
ond moment of the end-to-end vector of an ideal (non- 
interacting) a-type chain under flow, with its center-of- 
mass being in the ith layer. The last two terms in Eq. (JSJ 
are contributions from the kinetic energy and the free en- 
ergy change due to stretching of the chains. 

The conformations of interacting polymer chains are 
generated by lattice random walks in a self-consistent 
field modeling these interactions. According to the 
segmental volume fraction are obtained from the free seg- 
ment probability by 
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where 
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for s = 2 L 3, • • • , N a and Pf (1) = P?/(l - Pf - Pf). In 
Eq. is the average volume fraction of species a in 

the system and in Eq. Q , \± i is the stepping probability 
of a a-type segments in an ideal chain from layer i to layer 
i ± 1, whereas Aq i is the stepping probability from a site 
in layer i to a nearest neighbor in the same layer. In 
contrast to the static SCF theory, these probabilities are 
now time-dependent and anisotropic, and are related to 
the components of Sf through the random walk picture 
(see 0). 

Finally, Sf evolves according to 
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which is the constitutive equation of the FENE-P dumb- 
bell model 0, but with the value of the spring con- 
stant chosen in such a way that S" reverts to the 
Hookean dumbbell value in equilibrium. In Eq. JHJ), 
T% bs =N a2 Q a 2 /24:kBT is the local dumbbell relaxation 
time where £f is the local segmental friction coefficient of 
species a, which is related to the one in the bulk through 
the Doolittles' law (see for details). 

In the actual computation, all spatial derivatives are 
approximated by using finite difference formulas. Under 
suitable boundary conditions, Eqs. QJ, @ and © form 
a closed system of ordinary differential equations, which 
is solved numerically. 
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Henceforth we focus on symmetric blends, i.e. N A = 
N B = N, ni A = m B = m, and the two polymers A and B 
have the same density and segmental friction coefficient. 
We also set xf — xf — 0- We shall compare the numerical 
solutions of the DSCF equations with Molecular Dynam- 
ics (MD) simulations of a similar system. The details of 
the model and the method of the MD simulation that 
we used can be found in |8j. In MD, the characteristic 
energy (e), time (r) and length (c), and the segmental 
mass (m) are related by r = ayrnje. As in Q, we set 
the temperature in our MD simulations by e = fcsT/1.1. 

In order to compare our model to MD simulations, 
we have to establish the relationships between the model 
parameters in the DSCF theory and those in MD. The 
Kuhn length in the MD simulations is estimated to be 
a ~ 1.3cr 0, hence, one Kuhn segment is equivalent 
to 1.3 beads in MD. The interaction parameter xab is 
obtained by fitting the equilibrium density profiles cal- 
culated from DSCF theory to the ones obtained from 
MD. After accounting for capillary wave broadening of 
the intrinsic profiles in MD simulations [jj, this gives 
YIxabI^bT ~ 1.1. The segmental friction coefficient in 
the bulk is obtained by comparing the self-diffusion con- 
stants calculated from MD to the one in our model. The 
segmental mass in our model is identical to the one used 
in MD. In the simulations, we set (j> A = (f) B and the over- 
all polymer volume fraction averaged over the simulation 
box to (f> = 4> A + <J> B = 0.85. This average volume fraction 
is used both in the DSCF model and MD simulations. 

As the first test of the DSCF model, we calculated 
the steady state velocity profiles and the local viscosities 
between the walls under shear for different chain lengths, 
and compared them with the MD results. The steady 
state is obtained by integrating the DSCF equations of 
motion in time under a constant rate for shearing the two 
walls (nominal shear rate), starting from equilibrium. 

The shear velocity profile (defined as the local mass- 
averaged velocity) across the polymer-polymer interface 
in the middle of the channel from our model for N—12 
is shown in the left panel of Fig. ^ The result from MD 
simulations using chains of 16 beads is also shown. The 
shear rates 7 in the bulk are indicated in the caption |10| . 
This bulk shear rate is obtained by fitting the linear part 
of the velocity profile to a straight line and is different 
from the nominal shear rate because of the velocity slip 
at the walls and the interface. Note that a kink repre- 
senting a velocity slip at the interface is observed here. 
This interfacial slip phenomenon was first predicted by 
de Gennes |ll| based on a scaling argument, and later 
modeled with an approximate constitutive equation for 
incompressible polymer blends by Goveas and Fredrick- 
son IIS, and was subsequently confirmed by Barsky and 
Robbins 01 using MD. Experimental studies have also 
indicated the presence of a velocity slip at the interface 
in polymer blends 0|. Here the same velocity slip is 
produced in our DSCF model. The bulk shear rates for 
the DSCF and MD results shown in Fig. 2] are not iden- 
tical, but the closest that we can attain. This is because 
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FIG. 1: Left: Steady state shear velocity profile for 
N=12 obtained from DSCF at 7=2. 36x 10" 3 r _1 and MD at 
7=4. 04x 10 -3 t -1 . (The velocity is normalized by the veloc- 
ity at z = 10a.) Right: Variation of slip length with shear 
rate and with xab for N=12 obtained from DSCF. (Inset 
shows the slip length in a Vs XAB/ksT at nominal shear rate 
2.75x1CT 3 t- 1 .) 



_ (a) 








(b) 




















f 
















_ 












i 







f 


• 






j j 


3BEEEEQ..X. 

m 

Ea 


H-.-Q 

□ 


L 


0. 

/ •'^•!v:-H-;. 


b 



-20 -10 10 -20 -10 10 20 



z/o z/a 
FIG. 2: Second moment of the end-to-end distance for N—12 
obtained from (a) DSCF at 7=2. 36x 10 _3 r _ \ and (b) MD at 
7=4.04xl0~ 3 r- 1 . (+: S xx , x: S yy , O: S zz , □ : S xz .) 



simulating systems with low shear rates in MD is not 
computationally feasible. On the other hand, the DSCF 
simulations currently become unstable at higher shear 
rates, though they cover most experimentally accessible 
shear rates. 

To quantify the interfacial slip, we calculated the slip 
length, which is defined as twice the magnitude of the 
x-intercept of the straight line obtained from a linear 
fit to the linear part of the steady state velocity pro- 
file. The results are shown in the right panel of Fig. ^ 
Our model shows that the slip length decreases as the 
shear rate increases. We verified this result by perform- 
ing MD simulations of the same system at two different 
bulk shear rates: one at 4.04xl0 _3 r _1 and the other 
one at 8.27x1(R 3 t -1 . The slip lengths were found to 
be 1.95cr and 1.22er, respectively. The trend is the same 
as the one predicted by the DSCF model. We also cal- 
culated the slip length as a function of xab at a fixed 
nominal shear rate and the result is shown in the inset of 
the same figure. It shows that the slip length increases 
as the strength of the segment-segment interaction in- 
creases. A similar qualitative trend was also observed by 
using MD in earlier studies 0. 

We also calculated the second moment of the end-to- 
end vector at steady state by performing a Monte-Carlo 
simulation of lattice random walks in a self-consistent 
field representing conformations of interacting chains un- 
der shear. The transition rates used in the Monte-Carlo 
simulation were constructed from the stepping and free 
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FIG. 3: Steady state viscosity profiles across the polymer- 
polymer interface for different chain lengths obtained from (a) 
DSCF at nominal shear rate 6.24x 10~ 4 t~\ and (b) MD at 
various nominal shear rates. (+: N=6, x: iV=12, O: iV=18, 
□ : iV=24.) 

segment probabilities from the DSCF theory. The result 
is shown in Fig.[2]and compared to the one obtained from 
MD. In the bulk, we have good agreement between the 
two simulations. The slight difference can be accounted 
for by the difference in the shear rates. In the interfacial 
region, our model predicts the same trend as in MD. Sim- 
ilar results were also obtained by Barsky and Robbins in 
[3 using MD. 

The steady state viscosity profiles in the interfacial re- 
gion for various chain lengths are calculated by dividing 
the viscous shear stress by the local shear rate, and are 
shown in Fig.^a). Note that a drop in the viscosity oc- 
curs at the polymer-polymer interface and it is more pro- 
nounced for longer chains. Corresponding results from 
MD are shown in Fig. E^b). Here one can see that the 
trend predicted by the DSCF model is in general agree- 
ment with the MD results. The shear viscosities obtained 
from both methods in the bulk regions are very close. In 



the interfacial region, the drop is underestimated for long 
chains and slightly overestimated for short chains in the 
DSCF model, but the general agreement is very good 
given all the approximations of the model. A similar 
drop in the viscosity was also observed by Barsky and 
Robbins using MD 13], but with different parameters. 

We have studied the interfacial slip phenomenon in 
phase separated sheared polymer blends using the novel 
DSCF theory and compared our results with MD simu- 
lations. Good agreement between the two approaches is 
found. For unentangled polymer fluids studied here, the 
DSCF simulation time is only a small fraction of that 
in MD. We expect that these computational efficiences 
will be more manifest in the case of entangled polymer 
fluids, for which MD simulations are currently impracti- 
cal. The current DSCF theory is based on conservation 
laws for species occupancies and momentum that are cou- 
pled to models of polymer structure and conformation. 
It can be generalized either by changing the constitutive 
stress equation (accounting for the reptation regime), the 
polymer model (e.g., accounting for block copolymers or 
branching) or by incorporating additional conservation 
laws (e.g., for energy and/or charge). Work along these 
lines will be reported in future publications. 
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